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Abstract 

We study the interaction of an atom with a one-dimensional nano¬ 
grating within the framework of macroscopic QED, with special em¬ 
phasis on possible anisotropic contributions. To this end, we first derive 
the scattering Green’s tensor of the grating by means of a Rayleigh ex¬ 
pansion and discuss its symmetry properties and asymptotes. We then 
determine the Casimir-Polder potential of an atom with the grating. 
In particular, we find that strong anisotropy can lead to a repulsive 
Casimir-Polder potential in the normal direction. 

Keywords'. Casimir-Polder potential; Material grating; Rayleigh scat¬ 
tering. 
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1 Introduction 

Dispersion forces such as the Casimir-Polder interaction between micro¬ 
scopic particles (atoms and molecules) with macroscopic bodies are proto¬ 
typical forces that result from quantum vacuum fluctuations of a physical 
quantity, in this case the electromagnetic field mm- Their magnitude de¬ 
pends sensitively on the geometry of the macroscopic body and the optical 
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response of both particle and medium. In the simplest case of flat surfaces 
and isotropic optical response the Casimir-Polder force is usually attractive 
and perpendicular to the surface. However, in more complex geometries this 
force can be modified to a large extent. For example, surfaces with gratings 
can induce force components that are parallel to the (spatially averaged) 
interface. This has been shown for Casimir forces between corrugated sur¬ 
faces [3HS] as well as for Casimir-Polder forces [ZHin]. 

Moreover, the sign of the dispersion force can be reversed. A repulsive 
force with associated equilibrium, at least in the form of a saddle point, has 
been predicted for objects above a metal plate with a circular hole 
and near wedges m- Apart from geometry, non-equilibrium situations such 
as resonant Casimir-Polder interactions of excited atoms can provide tran¬ 
sient repulsion [MKIT!. In this article, we will show how a strong anisotropic 
optical response of an atom can lead to a repulsive Casimir-Polder force 
normal to a one-dimensional grating structure. 

From a theoretical point of view, to compute an electromagnetic disper¬ 
sion force means to determine the scattering properties of a surface, typically 
in terms of reflection matrices that enter a mode expansion mm or, equiv¬ 
alently, a Green’s tensor construction m- Here we use the dyadic Green’s 
tensor expansion in terms of Rayleigh reflection coefficients to investigate the 
interaction of an anisotropic atom with a one-dimensional periodic surface. 

2 Green’s tensor of a surface with one-dimensional 
periodic profile 

We require the scattering Green’s tensor (r, r', a;) for source point r' 
and field point r being situated in the vacuum half space above a nano¬ 
grating {y,y' > 0) (see Fig. [1]). The grating displays a periodic surface pro¬ 
file f{x) G [—/i, 0] in the x-direction with period d, f{x + d) = f{x) and 
height h while being translationally invariant in the z-direction. To con¬ 
struct the Green’s tensor, we employ a plane-wave basis which is adapted 
to the symmetry of the system. Wave vectors for upward (-I-) and down¬ 
ward (-) moving waves are parametrised as k± = (/c™, Their x- 

components are decomposed according to = kx + rnq into a continuous 
part kx G [—q/2,q/2\ that lies within the first Brillouin zone and a discrete 
part comprised of integer multiples m G Z of the reciprocal lattice vector 
q = 21: jd. The z-component takes arbitrary continuous values G M while 
the y-component then follows from the dispersion relation: 

/c™ = k'^{kx,kz,uj) = lim ^ (1 -|-i(5) — (A:™)^ — k"^ with ImA;™ > 0. 

( 1 ) 
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Figure 1: An atom (black ellipse) at position (0, y, z) above a one¬ 

dimensional rectangular grating. 


For each such wave vector k± = km±{kx, kz,oj), we choose two perpendicular 
polarisation vectors and such that only the respective electric and 
magnetic helds exhibit non-vanishing y-components, respectively: 


e^^{kx,kz,u) = 


_Um-L 

nih'^h 


^m±{kx: kz, to) — 


co^ioyc'^ - kl ’ 
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( 2 ) 

( 3 ) 


These vectors are obviously normalised, = 1 and or¬ 
thogonal to each other, = 0, as well as to the unit wave vector 

gfcj- = k±/\k±\, = 0. The three orthonormal vectors 

form a right-handed triad, e^j_ xe®j_ = efc±. 

With these preparations at hand, the required scattering Green’s tensor 
can be given as 


■ l^q/2 r 

G(‘)(r,r',i.) = ^/ dfe 5: I 

9/^ m,n=-oo"' 


dA:, 


E 


a,a'=E,H 

^i[k^x-k^x'+kl^y+k^y'+kz{z-z')] 
Un 

Ky 


'^m+ ^ '‘mn '^n— • 


( 4 ) 


The Rayleigh reflection coefficients — Rmnikx,kz,to) conserve kx and 
kz of the incident wave, but they mix polarisations a, a' as well as diffraction 
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orders m,n |2n[l21| . In general, they have to be calculated numerically 
by integrating the Maxwell equations within the grating {—h < y < 0) and 
imposing conditions of continuity at its upper and lower boundaries [3H5|I^. 
For a rectangular grating (Fig. [T|), the formalism for computing the Rayleigh 
coefficients is presented in[5l Before we proceed, let us derive some symmetry 
properties of the reflection coefficients as well as consider the asymptotes of 
small and large grating periods. 


2.1 Symmetry properties 

Schwarz reflection principle. Like every causal response function, the scat¬ 
tering Green’s tensor obeys the Schwarz reflection principle —oj*) 

= G*'^^*(r,r',w) [1]. Applying this to Eq. ([3]) leads to 

= R^r^nikx^h^io). (5) 


Onsager reciprocity. Provided that the grating consists of a medium 
which obeys time-reversal symmetry [22], the scattering Green’s tensor fulfils 
Onsager reciprocity G^^^(r',r,a;) = r',a;). This implies 


kx, kz,U}) 


Um 

Ky 


^ ^ Rml{kx,k^,uj) 
kV: 


for a = a',a a'. (6) 


Grating symmetries. Finally, we consider symmetries imposed by the 
geometry of the grating. The invariance of the grating with respect to an 
inversion of the 2 : coordinate implies that the diagonal/off-diagonal Rayleigh 
coefficients are even/odd functions of kz, respectively: 

Rmnikx,-kz,^^) = ±Rmnikx,kz,uj) for (T = cr',fj / ct'. (7) 


For a symmetric grating with /(—x) = /(x), we further have 


R’flCi-ni-kx,kz,oj) = ±R'^^{kx,kz,uj) for a = a',a a' 


( 8 ) 


by virtue of the x-inversion symmetry. 

Further symmetries can be obtained by combining the above. For in¬ 
stance, Eqs. dSI) and (l6|) imply R'^^*{kx,kzfl^)/k'^ = ±R^^{kx,kzflO/ky 
for a = a',a ^ a' at purely imaginary frequencies cj = i^. 


2.2 Asymptotes 

The exponential factors ^ and become either exponentially damped 
or rapidly oscillating whenever \k'^\y > 1 or \ky\y' > 1, respectively. In the 
limit where the grating period is small with respect to the distances of source 
and field points from the grating, d y,y', this implies that the Rayleigh 
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sums in Eq. 0 are effectively limited to small values. To leading order 
m = n = 0, we have 


^ik°{y+y')+ik^{z-z') 




in the coincidence limit r = r'. This asymptote of the scattering Green’s 
tensor is thus translationally invariant along the x-direction and it is de¬ 
scribed by the effective properties of the grating averaged over one period. 
The leading correction to this x-independent Green’s tensor is a harmonic 
variation in qx which is due to terms m = 0, n = ±1 and n = 0, m = ±1. In 
particular for a symmetric grating, this correction is proportional to cos{qx) 
in the coincidence limit. 

In the opposite extreme of the grating period being very large with re¬ 
spect to the distances of source and field points from the grating, d S> y, y', 
very large Rayleigh reflection orders m, n 3> 1 contribute. We thus have 

^ mq, k'^{kx,kz,oj) Ri k'^{0,kz,uj), so that the polarisation unit vec¬ 
tors and reflection coefficients can be approximated by e'^j.{kx,kz,oj) ^ 
e5^^(0, kz,uj) and R'^^{kx, kz,uj) ~ R^^{0, kz,uj)- Garrying out the remain¬ 
ing A: 3 ;-integral, the scattering Green’s tensor assumes the asymptotic form 



( 10 ) 


In the coincidence limit, the scattering Green’s tensor is then well described 
by a proximity force approximation |23ll24j where the grating is replaced by 
the local surface at x = x^ 

3 Casimir—Polder potential above a surface with 

one-dimensional periodic profile 

Within leading-order perturbation theory, the CP potential of a ground- 
state atom within an arbitrary structure is given as [25] 



( 11 ) 
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where 


^/cO^O/c 


^Ok^kO 


( 12 ) 


a{uj) = lim — , 

e->o+ n V w + Wfc + le a; — Wfc + le 

k ^ 

is the ground-state polarisability of the atom as given in terms of its transi¬ 
tion frequencies ujk = {Ej^ — EQ)/h and dipole matrix elements dmn- Using 
the Green’s tensor (Uj) from the previous section, we find the CP potential 
of an atom above a one-dimensional grating as 


Uir) = 


h/ip 

IbTT^ 


de^2 


rq/2 


^iq{m—n)x 


r 




■^—{Km+K-n)y 


^ d/c^ e'*? 

m,n=—oo 

X (13) 

a,a'=E,H 


with = R^^{k^,kz,iO and 


Km = Km{kx, k^, i^ = \l K + (k^f + • 


(14) 


In the limit d <C y of small grating period, we use the asymptote dH) of 
the Green’s tensor to hnd the x-independent potential 




dkz 




f —oo ^0 

X Y ^oo'eo+-Q:(iO-eol. (15) 

cr,a'=E,H 


For large grating period d ^ y, Eq. (flU]) leads to 


Uir) = 


kjiQ 

Svr^d 


poo pc 

/ Y / 

0 -TV. —_rv^ J 




m,n=—oo 


dk.- -- y; R^mi 

a,a'=E,H 

X e^+-Q(iO-e^'_ (16) 

with kx = 0 taken in all expressions inside the integrand. 


4 Casimir—Polder potential above a one-dimensional 

rectangular grating 

We study a ground-state Rb atom above a rectangular Au grating as shown 
in Fig. [TJ The grating consists of a periodic array of rectangular bars of 
height h = 20 nm and width 2 ym with separation di = 2 ym, so that the 
grating period is d = 4 /xm. The atom-grating potential is derived by nu¬ 
merically integrating Eq. ()13p with the Rayleigh coefficients as given in[5j 
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Figure 2: Contributions to the Casimir-Polder potential of a ground-state 
Rb atom above an Au grating as shown in Fig. [TJ 

To study the effect of possible anisotropies on the potential, we separate 
the latter into contributions from axx{'^) = oiyy{uj) = ayy{oj) = a{ui) where 
a{uj) is the ground-state polarisability of the Rb atom: U{r) = Uxx{f')+ 
Uyy{r) + Uyy{r). Tho resulting potentials at x = 0 are shown in Fig. [2j We 
find that the component Uyy leads to a repulsive Casimir-Polder potential in 
normal direction at separations y < 700 nm. This is the first demonstration 
of a repulsive Casimir-Polder potential for a grating geometry. 

In many cases of interest, a full numerical calculation of the potential 
is neither desirable nor needed. For short atom-grating separations, the 
asymptotic expression (jlOp gives reliable results. In Fig. [3l we show a com¬ 
parison of the full numerical calculation according to Eq. (I13|) (thick lines) 
with the short-distance asymptote (I16p (thin lines) at a fixed atom-grating 
distance of y = 700 nm. The obvious appearance of a non-sinusoidal lateral 
force across the grating implies in particular that the repulsive Uyy compo¬ 
nent at X = 0 leads to an unstable, saddle-point equilibrium in accordance 
with the generalised Earnshaw theorem [26]. One further notices that the 
short-distance (or large grating) approximation already provides very accu¬ 
rate results already at y = 700 nm. At smaller separations of y = 300 nm, the 
difference between exact and approximate potentials are within the thickness 
of the lines. 
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Figure 3: Casimir-Polder potential of a ground-state Rb atom at fixed 
distance y = 700 nm above an Au grating as a function of lateral position. 
Thick lines: exact numerical results from Eq. (m. Thin lines: asymptotes 
from Eq. (fTOt) . 

5 Conclusions 

We have shown that it is possible to generate repulsive Casimir-Polder forces 
above material gratings using anisotropic atoms. Our calculation was based 
on the Green’s tensor expansion of the electromagnetic field which, using a 
Rayleigh expansion, yields both a numerically exact result as well as ana¬ 
lytically tractable asymptotes for small and large gratings. Surprisingly, the 
short-distance (or large grating) asymptote becomes already accurate at a 
distance of 700 nm above a grating with a 4 ym. period. 
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Appendix A. Rayleigh reflection coefficients for a 
rectangular grating 

In the following, we determine the Rayleigh reflection coefficients for the 
rectangular grating depicted in Fig. [TJ To this end, we express the electro¬ 
magnetic fields above and below the grating for incident waves of polarisa¬ 
tion o' according to Eqs. Q and ([3]) as (the factor e is omitted for 

brevity) 


Ez{x,y) 

CO 

(17) 

Hz(x,y) 

n=—oo 

oo 

(18) 

for y > 0 and 

n=—oo 

oo 

E,{x,y) = Y 

(19) 


n=—oo 

oo 

H,(x,y) = Y 

n=—oo 

(20) 


for y < —h, respectively. For simplicity, we take c = 1 throughout this Ap¬ 
pendix. 

In the region —h < y < 0 one can write Ez{x,y) = Yl'^=-oo iv) 
exp(ifc”x), analogous decompositions hold for the other components of the 
electromagnetic field. It is convenient to denote by [E] the 2AI-|-1-component 
vector with components E]\f, Ej\f-i, ..., Eq, ..., E_jv+i, E_Ar. We further in¬ 
troduce a diagonal {2N -|- 1) x {2N + l)-matrix 

A = diag(/c^, /cf ..., ..., (21) 

The Maxwell equations in the region —h<y<0 can then be written as |21j 

® - ik,[Ey] = iu;[H^], ® - ik,[Hy] = -iu[D^], (22) 

ik:,[Ea;] - iA[Ez] = iuj[Hy], i/c^[i?a;] - = -iuj[Dy], (23) 

iA[Ey] - ^ = iu:[H,], iA[Hy] - ^ = (24) 

Writing [D] = Q[E] and using the results of Ref. [27|, one has for a rectan¬ 
gular grating: 
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Here, ||e|| is a Toeplitz matrix defined as 


/ eo 

£l 

£2 

^2N \ 

£-1 

£0 

El 

£2Af-l 

£-2 

£-1 

Eq 

£2N-2 

\£-2Af 

£-2iV+l 

£-2N+2 ■■ ■ 

£o / 


where En is a Fourier coefficient of the periodic function e{x + d) = e{x): 

^ exp(^^^^^^e„. (27) 


In order to solve the above system of equations, we proceed as follows. 
Equation (p3]l implies [Ey] = ||e||“^(A[77;j] — kz[Hx\)/ijJ., which can be sub¬ 
stituted into equations (l2^ and to eliminate [Ey], Similarly, Eq. (123)1 
leads to [Hy] = {kz[Ex] — A[Ez])/uj, which can be substituted into equations 
(j22p and (j24ll to eliminate [Hy], As a result, the Maxwell equations in the 
region —h<y<0 can be rewritten as 


/[A]\ 


/[E.]\ 

[A] 

= M 

[Ez] 

[H.] 


[H.] 



\[Hz]/ 


(28) 


where the matrix M is defined as 


M = i 


0 

0 


^A 


\ k 1 T I 


UJ 


1 ll-l 


0 

0 

kll - 

-%A 


-^fAikir' 

-%\\e\\-^+u:I 

0 

0 



-^A-cjA 

e||-A 


0 

0 


/ 

(29) 


Knowing the transmitted fields (fTHD and (I^UD at y = —h, one can deter¬ 
mine the fields at y = 0 by integrating Eq. (I28p from —h to 0. Imposing 
the continuity conditions on all Fourier components of [E^], [Ez], [H^], [Hz] 
at y = 0, one then determines the Rayleigh reflection and transmission 
coefficients. 
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